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1 Introduction 

Quantum Chromodynamics (QCD) is well established as theory of strong 
interactions through a large number of experimental verifications. The era 
of ‘testing QCD’ is clearly finished, and QCD today is becoming precision 
physics. In contrast to LEP and SLC, where electroweak reactions could 
be studied to a highly precise level without having to take QCD effects into 
account, QCD is ubiquitous at hadron colliders, affecting all observables. Any 
precision measurement (strong coupling constant, quark masses, electroweak 
parameters, parton distributions) at the Tevatron and the LHC, as well as 
any prediction of new physics effects and their backgrounds, relies on the 
understanding of QCD effects on the observable under consideration. 

Precision QCD poses several severe challenges from the theoretical and 
computational point of view. Most importantly, the strong coupling constant 
is considerably larger than the electromagnetic coupling constant at scales 
typically probed at colliders: as{Mz) — 15aemiMz), resulting in a slower 
convergence of the perturbative expansion. As a consequence, a precise de¬ 
scription of QCD observables (precise means here that the theoretical un¬ 
certainty becomes similar to the achieved or projected experimental errors) 
is obtained only by including higher order corrections, often requiring be¬ 
yond the next-to-leading order. Another important challenge is the fact that 
QCD describes quarks and gluons, while experiments observe hadrons. This 
mismatch is either accounted for by a description of the parton to hadron 
transition through fragmentation functions or by defining sufficiently inclu¬ 
sive final state observables, such as jets. Finally, many collider observables 
involve largely different scales, such as quark masses, transverse momenta 
and vector boson masses. These give rise to potentially large logarithms, 
which might spoil the convergence of the perturbative series and need to be 
resummed to all orders. 

In this talk, I shall try to highlight recent theoretical progress towards pre¬ 
cision QCD at colliders, focusing on heavy quark production in Section H, on 
jets and multiparton final states in Section ^ and on photons and electroweak 
bosons in Section Finally, a summary of the current state-of-the-art and 
of yet open issues is given in Section p. 

* Invited talk at the 14 th Topical Conference on Hadron Collider Physics 
(HCP2002), Karlsruhe, Germany, 30 Sep-4 Oct 2002 
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2 Heavy Quarks 

The production of heavy quarks is one of the main topics investigated at 
high energy collider experiments. Heavy quarks are of particular interest to 
elucidate the flavour sector of the standard model, which is less well tested 
than the gauge sector. Also, many approaches to physics beyond the standard 
model, often related to electroweak symmetry breaking and mass generation, 
predict new effects to be most pronounced in observables involving heavy 
quarks. 

2.1 Total Cross Sections 

The total cross sections for the production of heavy quarks can be computed 
within perturbation theory. The current state-of-the-art is a next-to-leading 
order (NLO) calculation |^, which is further improved by summing large 
logarithms due to soft gluon emission up to the next-to-leading logarithmic 
(NLL) level ||]. As can be seen from Figure |^, these predictions are in good 
agreement with experimental data on the total tt cross section at the Teva- 
tron [H and the total bb cross section at HERA-B (which both actually 
refer to similar kinematical values of mq / y/s). The theoretical uncertainty on 
the prediction for HERA-B is larger for two reasons: the larger value of the 
strong coupling at mt than at to* and the dominance of gg initial states in pN 
collisions (HERA-B) compared to qq dominance in pp collisions (Tevatron). 

2.2 Transverse Momentum Distributions 

Differential distributions of hadrons containing b quarks measured in hadron- 
hadron, photon-hadron or photon-photon collisions have been in apparent 
discrepancy with theoretical predictions for quite some time. The spectrum 
of hadrons measured at CDF [|| is one of the most recent examples for 
this discrepancy. 



Fig. 1. Total cross sections for tt at the Tevatron and bb at HERA-B. 
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The theoretical prediction for B meson production involves a convolution 
of the hard matrix element for heavy quark production in parton-parton scat¬ 
tering with initial parton distributions and final state fragmentation functions 
describing the non-perturbative transition from a b quark to a i? hadron. It 
is in particular the latter which may account for the discrepancy between 
theoretical prediction and experimental data, especially since it has been 
observed that the transverse momentum distribution of 6-tagged jets [Q 
(which has little sensitivity to fragmentation functions) is in much better 
agreement with theoretical predictions. 

The definition of heavy quark fragmentation functions is not free from 
ambiguities, since some aspects of these functions are actually calculable in 
perturbation theory H]. In extracting these fragmentation functions from 
data on B hadron production in e+e“ collisions, several choices are made, 
related to the order of perturbation theory, the incorporation of mass effects 
in the matrix elements, the resummation of potentially large perturbative 
terms, the correction of data for parton showers or the parametric form of 
the ansatz used in the determination. A commonly used parametric form is 
the so-called Peterson ansatz [^, which contains only a single parameter e, 
describing the shape of the non-perturbative fragmentation function. Unfor¬ 
tunately, the sensitivity of the fragmentation function on the assumptions 
used in the extraction from e+e“ spectra is often overlooked when using this 
fragmentation function to compute heavy hadron spectra at colliders. It was 
pointed out in 0 that a consistent treatment (in this case a purely mass¬ 
less approach) in extracting and implementing fragmentation functions can 
yield a sizable reduction of the discrepancy between experimental data and 
theoretical prediction. 

In the kinematical range covered by the heavy hadron production at 
Tevatron, b quark mass effects are sizable II . An approach incorporating 
quark mass effects, perturbatively calculable components of the heavy quark 
fragmentation function and resummation of large logarithms up to the 
next-to-leading logarithmic level was presented in [|ll| with the fixed-order 
next-to-leading log (FONLL) scheme, which requires only a small, genuinely 
non-perturbative component of the fragmentation function to be fitted to 
e^e“ data. In order to expose the information content actually relevant to 
heavy hadron spectra at hadron colliders, this fit is done in moment space. 

In view of new data from ALEPH |^, a phenomenological study of B 
hadron production at colliders based on the FONLL scheme was pre¬ 
sented in . It was shown that the consistent treatment of the fragmenta¬ 
tion function in extraction and prediction reduced the discrepancy between 
data and theoretical prediction considerably, Figure ||. The theoretical predic¬ 
tion is however still falling somewhat short of the experimental data, which 
is probably due to currently uncalculated corrections beyond NLO. 

Many collider experiments also report an excess in the b quark produc¬ 
tion spectra. In interpreting these data, it must always be kept in mind that 
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Fig. 2. Transverse momentum spectrum of B hadrons at CDF, compared to calcu¬ 
lations using Peterson and FONLL fragmentation functions, from 


it is not b quark but B hadron production which is observed in the exper¬ 
iment. Information on b quark production is only inferred from these data 
using some model for the heavy quark fragmentation. As discussed above, 
there are numerous ambiguities, which can yield inconsistent predictions if 
not implemented consistently. In view of the rather sizable effects due to a 
consistent treatment of the fragmentation function observed on the B hadron 
spectra at CDF, it might be that the data sets on b quark spectra have to be 
reanalysed incorporating the new experimental information on the b quark 
fragmentation functions in a consistent manner. 

2.3 Top Quark Spin Correlations 

Contrary to b and c quarks, which have lifetimes much longer than the typical 
timescale required for hadronization, t quarks are too short-lived to hadronize 
before they decay. As a consequence, top quark production and decay can 
be fully computed within perturbative QCD, and the spin state of the t 
quark determined from the distribution of its decay products reflects 
the t quark spin induced in the production process, giving rise to non-trivial 
spin correlations in ti pair production processes [T^ . At present, data from 
Tevatron Run I are not yet sufficiently precise to probe these spin correlations, 
results obtained up to now concern only the polarization state of the vector 
boson produced in the t decay . Future measurements at Run II and LHC 
will allow to use tt spin correlations as a sensitive probe of potential effects 
beyond the standard model showing up in the electroweak vertex of t decay. 

3 Jets 

Hadronic jets at large transverse momenta are produced very copiously at 
colliders. Final states with a small number of jets are measured to very high 
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Fig. 3. Single jet inclusive cross section at NLO compared to DO data 
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Fig. 4. Dependence of the NLO prediction for the 3? cross section at the Tevatron 
on renormalization and factorization scale, from |l 


experimental accuracy, such that they can be used for precision measure¬ 
ments of the strong coupling constant and of parton distribution functions. 
Multiparton final states, involving a large number of jets, can on the other 
hand mimic final state signatures introduced by physics beyond the standard 
model, thus forming an irreducible background to searches. 

The current state of the art for the theoretical description of jet observ¬ 
ables at hadron colliders is next-to-leading order QCD. To this order, parton 
level event generators are available for Ij + X (EKS Q and JETRAD fisll), 
2j + X (JETRA D ra ), Sj + X (NLOJET-h-f- fll) andforE-hX, V+lj + X 
(both DYRAD ^), V + 2j + X (MCFM |^. In particular, the MCFM 
project aims to provide a library of a large number of hnal states to NLO ac¬ 
curacy. Final states with more than three identified particles (vector bosons 
or jets) are at present only known to leading order. NLO predictions do in 
general yield a very good description of the experimental data, e.g. for the 
Ij inclusive cross section [Q, Figure ||. The uncertainty on the theoretical 
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prediction, as estimated by varying renormalization and factorization scales, 
is reduced considerably from LO to NLO, Figure |^. Nevertheless, for observ¬ 
ables such as the single jet inclusive cross section, the theoretical errors of 
the NLO prediction are still in excess of the experimental errors. 

3.1 Jet Definitions 

Jets in high energy physics experiments are identified by applying a jet algo¬ 
rithm to hadronic final states. The default jet algorithm used in the analysis 
of Tevatron Run I data (iterative fixed cone algorithm) turned out to be un¬ 
suited for the study of multi-jet cross sections and to be infrared unsafe if a 
theoretical implementation beyond NLO is attempted. These problems can 
essentially be truncated back to difficulties in splitting and merging neigh¬ 
bouring cones, which is done using an ad-hoc separation parameter. Several 
improved algorithms have been suggested, which can be grouped into two 
classes: the fcr type algorithm |2^ uses a clustering procedure combining 
neighbouring preclusters according the transverse momentum of each clus¬ 
ter and to their mutual spatial separation. The improved legacy cone algo¬ 
rithm is an improved version of the iterative fixed cone algorithm, which 
avoids the problems of the latter. 

3.2 Precision Jet Physics 

Despite the evidently good agreement of NLO QCD with experimental data 
on jet production rates. Figure ^ predictions to this order are insufficient for 
many applications. For example, if one uses data on the single jet inclusive 
cross section j2^ to determine the strong coupling constant as, it turns out 
that the dominant source of error on this extraction is due to unknown higher 
order corrections. Given that the theoretical prediction to infinite order in 
perturbation theory should be independent of the choice of renormalization 
and factorization scale, this error can be estimated from the variation of the 
extracted as under variation of these scales, as seen in Figure |[ As a result, 
CDF find from their Run I data 

as{Mz) = 0.1178 ± 0.0001(stat)t;(;!((J^5(sys) ;(j ;(!(^7(scale) ± 0.0059(pdf). 

It can be seen that the statistical error is already negligible; improvements 
in the systematic error can be anticipated in the near future. To lower the 
theoretical error, it is mandatory to compute next-to-next-to-leading order 
(NNLO) corrections to the single jet inclusive cross section. 

A similar picture is true in e^e“ annihilation into three jets and deep 
inelastic (2 -|- 1) jet production, where the error on the extraction of ag from 
experimentally measured jet shape observables is completely dominated 
by the theoretical uncertainty inherent in the NLO QCD calculations. 

Besides lowering the theoretical error, there is a number of other reasons 
to go beyond NLO in the description of jet observables [^. While jets at 
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Fig. 5. Errors on extraction of Os from single jet inclusive cross section at CDF . 


NLO are modelled theoretically by at most two partons, NNLO allows up 
to three partons in a single jet, thus improving the matching of experimen¬ 
tal and theoretical jet definitions and resolving the internal jet structure. At 
hadron colliders, NNLO does also account for double initial state radiation, 
thus providing a perturbative description for the transverse momentum of 
the hard final state. Finally, including jet data in a global NNLO fit of par- 
ton distribution functions, one anticipates a lower error on the prediction of 
benchmark processes at colliders. 

The calculation of jet observables at NNLO requires a number of different 
ingredients. To compute the corrections to an n-jet observable, one needs the 
two-loop n parton matrix elements, the one-loop n-\-l parton matrix elements 
and the tree level n + 2 parton matrix elements. Since the latter two con¬ 
tain infrared singularities due to one or two partons becoming theoretically 
unresolved (soft or collinear), one needs to find one- and two-particle sub¬ 
traction terms, which account for these singularities in the matrix elements, 
and are sufficiently simple to be integrated analytically over the unresolved 
phase space. One-particle subtraction at tree level is well understood from 
NLO calculations |2^ and general algorithms are available for one-particle 
subtraction at one loop Tree level two-particle subtraction terms have 
been studied in the literature , their integration over the unresolved phase 
space was up to now made only in one particular infrared subtraction scheme 
in the calculation of higher order corrections to the photon-plus-one-jet rate 
in e^e“ annihilation A general two-particle subtraction procedure is still 
lacking at the moment, although progress on this is anticipated in the near 
future. 

Concerning virtual two-loop corrections to jet-observables related to 2 ^ 
2 scattering and 1 —*■ 3 decay processes, enormous progress has been made in 
the past two years. Much of this progress is due to several technical develop- 
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Fig. 6. Computer algebra to compute two-loop scattering matrix elements for jet 
physics. 


ments concerning the evaluation of two-loop multi-leg integrals. Using itera¬ 
tive algorithms |]^ , one can reduce the large number of two-loop integrals by 
means of integration-by-parts and Lorentz invariance 134 identities to a 
small number of master integrals. The master integrals relevant to two-loop 
jet physics are two-loop four-point functions with all legs on-shell or one 
leg off-shell , which were computed using explicit integration or implicitly 
from their differential equations . 

Combing the reduction scheme with the master integrals, it is straight¬ 
forward to compute the two-loop matrix elements relevant to jet observables 
using computer algebra The generic structure of such a calculational 
procedure is depicted in Figure Following this procedure, two-loop matrix 
elements were obtained for Bhabha scattering , parton-parton scattering 
into two partons @ , parton-parton scattering into two photons @ , as well 
as light-by-light scattering . Most recently, two-loop corrections were com¬ 
puted for the off-shell process 7 * ^ , relevant to e+e“ —> 3j. Part of 

these results were already confirmed using an independent method 1 ^ . 
Related to e+e” ^ 3j by analytic continuation are {2 + l)j production in 
ep collisions and V + j production at hadron colliders. A strong check on all 
these two-loop results is provided by the agreement of the singularity struc¬ 
ture with predictions obtained from an infrared factorization formula |46|. 


3.3 Multiparton Processes 

Many scenarios for physics beyond the standard model predict visible signals 
only in observables with a large number of partons in the final state due to 
production and subsequent decay of yet unknown short-lived particles. The 
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same multi-parton final states can however also be produced through ordi¬ 
nary QCD processes, acting as a background to the new physics signals. The 
extraction of new physics effects from multiparton final states does therefore 
require predictions for QCD backgrounds for these final states. 

It turns out that a diagrammatic evaluation of the corresponding matrix 
elements (even using modern techniques such as helicity amplitudes) becomes 
rapidly infeasible due to the factorial increase in the number of Feynman 
diagrams with increasing number of final state particles. An alternative ap¬ 
proach circumvents the need for a direct evaluation of Feynman diagrams 
by determining scattering matrix elements directly from the numerical Legen¬ 
dre transformation of the interaction Lagrangian. This procedure yields only 
a power-like growth of computation time with increasing number of final state 
particles, and does therefore outpace the Feynman diagrammatic evaluation 
for high multiplicity final states. Extending earlier work [Q on final states 
with up to nine gluons, a parton level event generator ALPHGEN to 
compute 

• VF/Z/y* -I- QQ + {n< 4)j 

• H + QQ + {n < 4) j 

• WiZh* + [n< 6)j 

• nW + mZ IH + Nj{n + m + l + N<8,N <3) 

• QQ + (n ^ 6)j 

• QQ + Q'Q' + {n< 4)j 

was presented recently. The event generation provides full colour and flavour 
information, such that a hadronization model and a subsequent detector sim¬ 
ulation can be interfaced. 

Combining matrix element based event generation with parton showers, 
special care has to be taken to avoid overcounting of contributions. Extending 
a method derived for e“'"e“ annihilation 1 ^ , a procedure for reweighting of 
events and determining appropriate starting conditions and cut-off of the 
parton shower has been proposed ||5l| to allow combination with the exact 
multi-parton matrix elements. 

A closely related issue is the incorporation of next-to-leading order cor¬ 
rections in parton shower Monte Carlo programs. Recently, algorithms were 
devised for this task by several groups ||52|, first programs implementing these 
algorithms have already been released |53|. 

4 Photons and Massive Gauge Bosons 

Photons and gauge bosons provide very prominent final state signatures at 
colliders. Their study allows the precise determination of electroweak param¬ 
eters at hadron colliders, and their final state signatures are often background 
to searches, such as photon pair production to the Higgs search in the lower 
mass range. 
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4.1 Isolated Photons 

Photons produced in hadronic collisions arise essentially from two different 
sources: ‘direct’ or ‘prompt’ photon production via hard partonic processes 
such as qg (77 and qq gj or through the ‘fragmentation’ of a hadronic 
jet into a single photon carrying a large fraction of the jet energy. The former 
gives rise to perturbatively calculable short-distance contributions whereas 
the latter is primarily a long distance process which cannot be calculated per¬ 
turbatively and is described in terms of the quark-to-photon fragmentation 
function. In principle, this fragmentation contribution could be suppressed to 
a certain extent by imposing isolation cuts on the photon. Commonly used 
isolation cuts are defined by admitting only a maximum amount of hadronic 
energy in a cone of a given radius around the photon. An alternative proce¬ 
dure is the democratic clustering approach suggested in Q, which applies 
standard jet clustering algorithms to events with final state photons, treating 
the photon like any other hadron in the clustering procedure. Isolated pho¬ 
tons are then defined to be photons carrying more than some large, predefined 
amount of the jet energy. 

Both types of isolation criteria infrared safe, although the matching of 
experimental and theoretical implementations of these criteria is in general 
far from trivial. It was pointed out recently that cone-based isolation 
criteria fail for small cone sizes R (once Oslni?”^ ~ 1), since the isolated 
photon cross section exceeds the inclusive photon cross section. This problem 
can only be overcome by a resummation of the large logarithms induced by 
the cone size parameter. 

The contribution from photon fragmentation to isolated photon cross sec¬ 
tions at hadron colliders is sensitive (for both types of isolation criteria) on 
the photon fragmentation function at large momentum transfer, which has 
up to now been measured only by one of the LEP experiments . Based on 
the democratic clustering procedure, ALEPH extracted the quark-to-photon 
fragmentation function from the measured photon-plus-one-jet rate. It was 
observed that the resulting prediction for the variation of the isolated 
photon rate with the resolution parameter was in good agreement with the 
measurement, especially once NLO corrections |]^ were included. A related 
OPAL measurement of the photon fragmentation function from inclu¬ 
sive photon production in e“'"e“ is unfortunately not sufficiently sensitive on 
the behaviour at large momentum transfers Eurther information on the 
photon fragmentation function at large momentum transfer might be gained 
from yet unanalysed LEP data or from the study of photon-plus-jet final 
states in deep inelastic scattering at HERA |^, where first data are now 
becoming available Q. 

4.2 Photon Pairs 

The discovery of a light Higgs boson {uih ^ 140 GeV) at the LHC is based 
largely on the observation of the rare decay to two photons . To perform 
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Fig. 7. Contributions from qq, qg and gg subprocesses to 77 final states for different 
isolation criteria, from |^. 
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an accurate background subtraction for this observable, one requires a precise 
prediction for QCD reactions yielding di-photon final states. At first sight, 
the 0{a^g) process qq 77 yields the leading contribution. However, due 
to the large gluon luminosity at the LHC, both qg — s- 977 (©(Ug)) and 
gg — > 77 (©(Ug)) subprocesses yield contributions of comparable magnitude 
(Figure ^). The NLO corrections to the qq and qg subprocesses have been 
known for quite some time, these are implemented in the flexible parton 


level event generator DIPHOX 
derived for the gg subprocess 


I . Most recently, NLO corrections were also 
Since the lowest order contribution to this 


process is already mediated by a quark loop, this calculation contains some of 
the features appearing in jet physics only at NNLO, such as two-loop QCD 


amplitudes and unresolved limits of one-loop amplitudes (see Section B.2 
above). 

It must be kept in mind that the di-photon cross sections are highly sen¬ 
sitive on the isolation criteria applied to the photons. Figure R, with a sub- 
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stantial contribution arising from photon fragmentation at large momentum 
transfers Moreover, it is experimentally difficult to distinguish photons 
from highly energetic neutral pions which decay into a closely collimated 
photon pair, mimicking a single photon signature. The pion background in 
photon pair production has been studied to NLO in and implemented in 
DIPHOX, showing that in particular the Tr^y channel remains comparable to 
the 77 channel even for tight isolation criteria. Figure |^. 


4.3 Vector Boson and Higgs Production 

The production cross sections for and bosons at hadron colliders are 
well understood both experimentally and theoretically. At present, these cross 
sections are measured to an error of about 10% from Tevatron Run I, largely 
limited by statistics. A considerable reduction of the experimental error is 
anticipated from Run II and for the LHC. On the theory side, inclusive vector 
boson production has been computed to NNLO already more than ten 
years ago. Very recently, these results have been verified for the first time in 
an independent calculation |67[. T he uncertainty on the theoretical prediction 
has been assessed in detail in ||68| and found to be around 3% from variations 
of renormalization and factorization scale as well as from uncertainties on 
the parton distribution functions. Some improvements are anticipated once 
full NNLO determinations of the parton distribution functions become 
available. 

Given the good theoretical and experimental understanding of and 

boson production, it has been suggested to use these for a determination 
of the LHC luminosity |^. In practice, it turns out that it is not possible to 
measure the fully inclusive production cross sections, but only cross sections 
integrated over a restricted range in rapidity, which are only known to NLO 
at present 0 . For the production, which is observed only through the 
decay channel, it is moreover mandatory to compute the spatial distribution 
of the decay products, which is again only known to NLO |^. To match the 
anticipated experimental accuracy at the LHC and to render vector boson 
production a reliable luminosity monitor, it would be required to extend both 
these calculations to NNLO. 

Closely related to the inclusive vector boson production cross section is 
the inclusive Higgs boson production cross section in the infinite top mass 
limit. While the former is induced, at the leading order, by the partonic 
subprocess gg —!■ W^,Z^, the latter is gg —> H. The NNLO corrections to 
inclusive Higgs boson production have been derived, first in the soft/collinear 
approximation |Q; shortly thereafter, the full coefficient functions were ob¬ 
tained by expansion around the soft limit and fully analytically by 
extending the IBP/LI reduction method and the differential equation tech¬ 
nique (see Section 3.2) to compute double real emission contributions. It 
turned out that inclusion of NNLO corrections yields a sizable enhancement 
of the Higgs production cross section, and a reduction of the uncertainty due 
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to renormalization and factorization scale. These calculations were also ex¬ 
tended to pseudoscalar Higgs production Q. A more detailed discussion of 
recent results on Higgs physics at colliders can be found in . 

4.4 Transverse Momentum Distributions 

Transverse momentum distributions of vector bosons or any other colour- 
neutral final states are well described in perturbation theory only if the 
transverse momentum qx is of the same order as the invariant mass M of 
the colour-neutral final state. Otherwise, the convergence of the perturba¬ 
tive series is spoilt by large logarithms h\{qT/M). Reliable predictions can be 
obtained if these logarithms are resummed to all orders in perturbation the¬ 
ory [f^ . This resummation procedure exponentiates logarithms from three 
different sources: soft (A) and collinear {B) radiation associated to the hard 
interaction (expressed as Sudakov form factor) and collinear initial state ra¬ 
diation (C). The resummation coefficients A, B, C can be expanded in powers 
of the strong coupling constant; it turns out that the soft coefficient A only 
depends on the partonic initial state, while B, C also depend on the final 
state under consideration. 

The leading logarithmic (LL) corrections are obtained from soft exponen¬ 
tiation alone, while collinear contributions first enter at the next-to-leading 
logarithmic (NLL) level 0- Finally, collinear initial state radiation affects 
only higher (NNLL) order logarithms. At present, NNLL resummed correc¬ 
tions are known for electroweak vector boson production and Higgs pro¬ 
duction . Incorporation of these corrections improves the theoretical 
description of vector boson production at small transverse momenta consid¬ 
erably, and yields a stabilisation of the theoretical prediction under variations 
of renormalization and factorization scales. 

The purely soft (A) terms can be read off from the partonic splitting func¬ 
tions to the appropriate order [6^ or computed by considering only eikonal 
diagrams Q. 

By reordering the terms in the resummation formula, it is possible to 
extract the final-state independent terms in the B and C coefficients, and 
to group the genuinely process dependent terms into a hard coefficient func¬ 
tion H. This modified resummation formula allows to obtain resummed 
predictions for a large number of different final states (such as vector boson 
pairs) without having to rederive the full set of resummation coefficients. 

In the course of the computation of these universal resummation 
coefficients to NNLL accuracy, one encounters configurations similar to the 
ones appearing in double real radiation contributions to jet physics at NNLO. 

Besides the above developments on extending the existing resummation 
formalism towards higher accuracy, considerable progress has been made also 
to perform resummation in processes not characterised by a colour neutral 
system with a large invariant mass. In particular, a formalism to jointly resum 
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partonic threshold and recoil corrections has been developed and applied 
successfully to prompt photon and vector boson production [Q. 

5 Conclusions and Outlook 

QCD at present and future hadron colliders will be precision physics, very 
much like electroweak physics was precision physics in the LEP era. The study 
of many of the standard scattering reactions will allow a precise determination 
of the strong coupling constant, electroweak parameters, quark masses and 
parton distribution functions. In turn, this information translates in improved 
predictions for new physics signals and their backgrounds. 

To match the experimental accuracy reached for a number of QCD collider 
observables, the current theoretical predictions have to be improved in sev¬ 
eral aspects. Most importantly, NLO QCD is often not enough if confronted 
with high precision data (such as for jets and gauge bosons) or observables 
with slowly converging perturbative expansions (such as heavy quarks or the 
Higgs boson). Considerable progress has been made very recently in extend¬ 
ing QCD calculations to NNLO, first results on inclusive observables are now 
available, and calculations for jet observables are well advanced. Many ob¬ 
servables do moreover require the resummation of large logarithms spoiling 
the convergence of the perturbative series. A universal picture for these re¬ 
summations started to emerge recently, and new techniques have been devel¬ 
oped to compute the corresponding resummation coefficients. Fragmentation 
effects enter many observables with identified particles in the final state. In 
particular, a consistent treatment of heavy quark fragmentation effects can 
account for a large part of the observed discrepancy in B hadron spectra, 
and quark-to-photon as well as quark-to-pion fragmentation yield important 
contributions to photon pair final states forming an important background to 
Higgs searches. Much valuable information on these fragmentation functions 
is contained in data from LEP, and should be extracted (as long as this is 
still a feasible task) to improve predictions for collider observables. 
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